Low temperature ballistic spin transport in the S—l/2 antiferromagnetic Heisenberg 

chain compound SrCuC>2 



(N 

o 

(N 
< 



H. Maeter/'B A. A. Zvyagin, 1 - 2 H. Luetkens, 3 G. Pascua, 3 Z. Shermadini, 3 
R. Saint-Martin, 4 A. Revcolevschi, 4 C. Hess, 5 B. Biichner, 5 and H.-H. Klauss 1 

1 Institute for Solid State Physics, TU Dresden, D-01069 Dresden, Germany 
2 Institute for Low Temperature Physics and Engineering of the NAS of Ukraine, Kharkov, 61103, Ukraine 
3 Laboratory for Muon-Spin Spectroscopy, Paul Scherrer Institut, CH-5232 Villigen PSI, Switzerland 
4 Laboratoire de Physico-Chimie de L'Etat Solide, ICMMO, 
UMR 8182, Universite Paris-Sud, 91405 Orsay, France 
5 IFW-Dresden, Institute for Solid State Research, P.O. Box 270116, D-01171 Dresden, Germany 

We report zero and longitudinal magnetic field muon spin relaxation (pSR) measurements of the 
spin S— 1/2 antiferromagnetic Heisenberg chain material SrCu02- We find that in a weak applied 
magnetic field Bo the spin-lattice relaxation rate A follows a power law Acx-B^™ with n=0.9(3). This 
result is temperature independent for 5 K<T<300 K. Within conformal field theory and using the 
Miiller ansatz we conclude ballistic spin transport in SrCu02. 

PACS numbers: 75.10.Pq, 76.75.+i, 73. 23. Ad, 76.60.Es 
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The copper oxide based low-dimensional electronic sys- 
tems allow a detailed study of many cooperative quan- 
tum phenomena and concepts like Mott insulators, spin- 
charge separation, quantum phase transitions, and un- 
conventional superconductivity 0, 0]. Quasi 1-d ar- 
rangements of corner-sharing Cu02 squares are model 
compounds for the S=l/2 antiferromagnetic Heisenberg 
chain (AFHC). In these systems charge degrees of free- 
dom are quenched at low energies by strong Coulomb 
interaction. Spin degrees of freedom are governed by the 
Heisenberg Hamiltonian H=JY). 5^5,;+! where i indexes 
the spins along the chain and the exchange constant 
J controls the interaction strength between neighboring 
spins. A model compound for the isotropic 5=1/2 AFHC 
is SrCu02- It is regarded as an almost ideal 1-d system 
with Jss2100 K between neighboring spins of magnetic 
Cu 2+ ions |3j. In relation to J the residual interchain 
interaction that causes magnetic order at temperatures 
below ?s2 K [4j is very small. 

The ground state (GS) of the S=l/2 AFHC is a sin- 
glet (5=0) state with a continuum of 5=1/2 excitations, 
spinons. Recently, using numerical solutions of the Bethe 
ansatz equations it was shown that the ground state dy- 
namic correlation functions are determined mostly by 
the 2-spinon continuum 0, 0] . While its thermodynamic 
properties have been studied thoroughly, both theoret- 
ically and experimentally, the dynamic properties, such 
as the spin transport have been much less studied: For a 
strongly interacting system with a non equidistant (due 
to interactions) spectrum of eigenstates it is an arduous 
task. As a result theoretical studies are often contro- 
versially discussed and even diffusive spin transport has 
been predicted for some temperature and field regimes 

Recent heat transport experiments revealed an unex- 
pectedly large magnetic contribution to the total heat 
conductivity of 5=1/2 AFHC materials such as SrCu02, 



attributed to the magnetic excitations of the spin chains 
at low temperatures 111, LL2J. It is still unclear how 



this contribution can be understood microscopically. 
A deeper understanding of the dynamic properties of 
the prototype Hamiltonian for the 5=1/2 AFHC, the 
Heisenberg Hamiltonian H— J 5^5^+1 is essential for 
progress in this field. Its most important property is in- 
tegrability, i.e. the existence of infinitely many local con- 
stants of motion. In gene ral, integrability implies bal- 
listic (spin) transport [l3|. In particular, ballistic spin 
transport is predicted for finite magnetic fields and the 
ground state jio) . However, experimental proof for bal- 
listic spin transport in the 5=1/2 AFHC is still lacking. 
Several experimental studies reveal diffusive spin trans- 
port in different model compounds of the 5=1/2 AFHC 
[l4H16[. In this Letter we present experimental evidence 
for ballistic spin transport in SrCu02 at low temperature 
T^C J and finite fields i?<§C J based on /LtSR experiments. 

The structure of SrCu02 contains chains running along 
the crystallographic c axis build of corner sharing CuO 
squares reminiscent of the CuO layers of cuprate super- 
conductors. In SrCu02 two chains are joined by shar- 
ing their edges, forming a zig-zag chain. The antifer- 
romagnetic coupling between nearest neighbors (NN) is 
large compared to the ferromagnetic exchange coupling 
J'«220 K between diagonal Cu spins [17] . Hence, it can 
also be described as a 5=1/2 chain with ferromagnetic 
NN and antiferromagnetic next nearest neighbor (NNN) 
interactions. The resulting weak frustration can produce 
essential features in the behavior of the spin chain. In- 
elastic neutron scattering data (INS), however, indicates 
that both chains are decoupled and show no features as- 
sociated with frustration Neutron diffraction exper- 
iments revealed anisotropic spin freezing below «5 K in 
SrCu02 [Hj]. This frozen state is also detected by /LtSR, 
but only below ~2 K [1] . Thurber and coworkers studied 
the dynamics of the q=0 modes in SrCu02 by 17 nuclear 
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FIG. 1. Typical uSR time spectra of SrCu02 at a temperature 
of 75 K measured for different longitudinal magnetic fields. 
Solid lines are best fits to the data (see text). 



magnetic resonance (NMR) [15}. They find a magnetic 
field Bq dependence of the spin lattice relaxation rate 
I/TiOcBq 71 that is consistent with diffusive spin trans- 
port, i.e. with n=0.5 [15j . 

In this study we used a single crystal of SrCuC>2 grown 
by the traveling solvent floating zone technique from high 
purity (99.99%) precursors. The thermal transport prop- 
erties of this sample have been reported by N. Hlubck 
and coworkers in Ref. fTlj j . A single crystal was oriented 
and mounted with the crystallographic b axis along the 
muon beam direction. /iSR experiments were conducted 
using a 4 He flow cryostat at the GPS instrument of the 
Paul Scherrer Institut, Switzerland. In a /iSR experi- 
ment, nearly 100% spin polarized muons are implanted 
into the sample one at a time. In cuprate materials the 
positively charged /i + usually form a ~1 A long bond 
with a oxygen ion where they act as magnetic micro- 
probes [l9j. In a non-magnetic material the muon spins 
will dephase, i.e., the 100% initial spin-polarization P{t) 
will decay as a function of time t due to the random ori- 
entation of nuclear magnetic dipole fields. The nuclear 
dipole field distribution can be modeled by an isotropic 
Gaussian distribution with width A. This is well known 
and can be described by the so called Kubo-Toyabe func- 
tion G(t, A, Bq), which, in a longitudinal magnetic field 
Bq has been described by Hayano [23]. Typically, £?o~5 
to 10 mT will decouple the muon spin from the nuclear 
dipole field distribution. 

The quantum spin fluctuations of the S=l/2 AFHC 
cause an additional independent relaxation mechanism 
for the muon spin ensemble. Rapid fluctuations will 
cause an exponential relaxation of the muon spin polar- 
ization P(i)ocexp(— At) with the spin lattice relaxation 
rate A. The overall relaxation function is then the prod- 
uct P(t)=G(t, A, Bo) exp(-Xt). In our experiments, A is 
of the order of 0.01 /is -1 and we study its field and tem- 
perature dependence. No signature of muon diffusion has 



been found, hence A=0. 085(2) /is -1 is temperature and 
field independent. Typical //SR time spectra and best 
fits to the data are shown in Fig. [TJ 

To confirm that the measurement is not influenced by 
the muon as an impurity we compare our data with recent 
63 Cu NMR data [2l| in Fig. [3] It turns out that both the 
/iSR and NMR spin-lattice relaxation rates have nearly 
identical temperature dependencies. This indicates that 
the muon has little or no influence on the low energy spin 
excitations of SrCu02- 

For diffusive spin transport in one dimension A oc 

13, 



22j. Its divergence for B o =0 can 
be cut-off by 3-d diffusion/coupling or by anisotropy, 
e.g. dipolar intrachain coupling. A useful description 
of experimental data has to include this cut-off for low 
fields/frequencies [231 ]: 



B 0,5 is expected 



X(B Q ) = A 



1 + y/1 + {B /2B c y 
2(1 + {B Q /2B C ) 2 ) 



(1) 



Ao=l/ y / 2D\\D± is the constant value obtained for 
Bq<B c and 7 /1 i3 c =Z?x, with n=0.5. Two diffusion con- 
stants D|| and D± model fast on-chain and slow intra- 
chain diffusion, respectively. We chose this model to 
maintain comparability with the work of Pratt et al. 14 1. 



Its main purpose is to determine the power law expo- 
nent n. In experiment 14], sometimes n^0.5 is found 
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FIG. 2. The field dependence of the relaxation rate \(Bo) 
shows a clear power law behavior for Bo>_B c =0.23(3) mT and 
saturates for smaller fields. The solid line is a fit of Eq. ([lj 
to the data (B c was optimized simultaneously for all temper- 
atures). The large deviation from the expected field depen- 
dence for _Bo=10 mT is due to the limited time window of the 
experiment of «f2 us. The open symbol indicates A(Bo=0). 

and (TJJ) becomes an empirical model since the deriva- 
tion in Ref. [23| is only valid for chains for ?i=0.5 (for 
a three-dimensional system, where -Dj_^0, the equation 
is formally correct), i.e. in the limit D±^0, X(Bo) van- 
ishes or diverges for n>0.5 or n<0.5, respectively. We 
use Eq. (fTJ) to analyze our data. For Bq>B c , X(Bq) fol- 
lows the power law -Bq - ™, here B c is called the cut-off 
field. A typical best fit of Eq. (TTJ to the data is shown 
in Fig. [21 £? c =0.23(3) mT is found to be temperature 
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independent, n also shows no temperature dependence 
for temperatures between 5 and 300 K, as can be seen 
in Fig. [3] The average value and its standard devia- 
tion are n=0.9(3). This result is quantitatively differ- 
ent from the result expected for diffusive spin transport, 
i.e. n=0.5. A microscopic interpretation of this result 
will be given below. Ao ~ 0.03 /js _1 could be identified 
with the inverse of a diffusion constant as in Eq. (pj and 
it shows no anomaly at low temperatures. In fact, all 
parameters that characterize the spin transport in this 
system, i.e. Ao, B c , and n have no significant tempera- 
ture dependence for T < 150 K. This indicates that the 
mechanisms that cause the temperature dependence of 
the magnetic heat transport, which shows a maximum 
at low temperatures 11], are not relevant for the spin 
transport. This corroborates the conclusions drawn in 
Ref. [ll|, that heat transport of the 5=1/2 AFHC in 
SrCuC>2 is limited, i.e. rendered diffusive by extrinsic 
perturbations, e.g. phonons. 
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FIG. 3. Top: A comparison of the temperature dependence 
of Ao with the spin lattice relaxation rate from Ref. [2l( mea- 
sured by 63 Cu NMR. The nearly identical temperature de- 
pendencies indicate that the muon does not disturb the low 
energy spin excitation spectrum. Bottom: The power law ex- 
ponent n. The shaded area indicates the standard deviation 
of n. All error bars are error estimates from the fits. 

Next, we will turn to the microscopic interpretation of 
the experimental results. The relaxation rate for a local 
probe, i.e. a nuclear, or muon spin in an electron spin 
system can be expressed as [24j 



+ -F^( q )[S + -(q,u;) + S- + ( q ,u;)}), (2) 

where 7 e (Ar) are the gyromagnetic ratios of electron (nu- 
clear) spins (2||, F z (q) and F^(q) are components of 
the hyperfine form factors parallel and perpendicular 



to the external magnetic field, S zz (q,uj) 5 H (<?,w), and 
5 ^(q,ui) are the components of the tensor of the dy- 
namic structure factor (DSF) of electron spins, also par- 
allel and perpendicular to the external field. In general, 
the muon relaxation rate A differs from the nuclear mag- 
netic resonance rate 1/Ti only by the form factor. 

A detailed knowledge of the form factor is not required 
because generally no "filtering" of the DSF by the form 
factor is expected in cuprate materials. This is due to 
the low symmetry of the muon site which is usually at a 
distance of «1 A from an oxygen ion [19j . 

The experimentally found \(xS zz (xl/ B^ with 
n=0.9(3) is in clear disagreement with the model for 
spin diffusion (n=0.5). Note that n=0.5 disagrees also 



with the theoretical work [26j of Sachdev [27[. This 
leaves us with two questions: (1) Is the spin transport 
ballistic?, and (2) How can we understand the observed 
power law? In the following we will show that the 
measured power law with n«l can be understood by 
the spin excitation spectrum of the ground state of the 
5=1/2 AFHC. In the ground state, exact calculations 
show that the spin transport is ballistic (see e.g. Ref. f5~oj ] 
and references therein). Therefore, nwl proofs ballistic 
spin transport in the 5=1/2 AFHC material SrCuC>2 
because we find n=0.9(3) by experiment. One model 
for the ground state spin excitation spectrum has been 
given by Miiller et al. This is a well known model and 
often used to analyze inelastic neutron scattering data 
1, 28] In this model S zz oc1/uj in agreement with our 
experimental result. We will describe this model below. 

According to the conjecture of Miiller et al., the ground 
state DSF of the 5=1/2 AFHC is determined by the 2- 
spinon continuum In the absence of the magnetic 
field it has a lower ei (q)=ir J/2 1 sing | and upper bound 
e2(q)=TrJsmq/2. The contribution to the ground state 
DSF of the 2-spinon continuum of the 5=1/2 AFHC can 
be written as 



S?(q,u) = A 



Q(m - e 1 (g))e(e 2 (q) - uj) 



(3) 



where A 2 is a constant, and O(x) is the Heaviside step 
function (See Ref. [29] for the optimized value of A 2 ). 

In the absence of a magnetic field the DSF is isotropic, 
and at uj— >wo^C J two points of the reciprocal space, <7=0 
and q=ir contribute mostly to the relaxation rate ([2]) [6] . 
On the other hand, a non zero external magnetic field 
introduces an anisotropy for the components of the dy- 
namic correlation functions of the 5=1/2 AFHC, shift- 
ing the contributing points away from g=0, it. However, 
in our /iSR experiments the external magnetic field was 
much smaller than the exchange constant along the spin 
chain, and this shift is negligible. In /iSR experiments 
we have uj=loq—^^Bq. In the ballistic regime, according 
to Eq. ([3]) the relaxation rate in the ground state has to 
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be inverse proportional to the value of the external field 



A K FT' 

-DO 



(4) 



Notice that the value Ai is field-dependent [30(. At low 
fields, this can cause a weak logarithmic dependence of 
the relaxation rate in addition to the inverse proportion- 
ality 0. The power law (gj), XocB^ 11 with n=l valid for 
the ground state DSF is in agreement with the experi- 
mentally found power law with n=0.9(3). According to 
above arguments, this is evidence for ballistic spin trans- 
port in SrCuC>2- In the following we will discuss the 
well known effects of exchange anisotropy, temperature 
and next nearest neighbor (NNN) exchange interaction 
on the power law ([4|). 

Uniaxial anisotropy of the exchange interaction can al- 
ter the frequency dependence of S zz (q,cu) and hence the 
exponent in the power law Q: 



S zz (q,co)=A 



ejiJ - e 1 (g))e(e 2 (q) - lj) 
(u 2 -el(q)) a °(e1(q) - w a)(l/2)-a 



(5) 



The exponent ao can be calculated from the finite size 
corrections to the ground state energy and a>o=l — 
Z 2 =(vr - 2?7)/(27r - 2rj), where Z=y/n/2(n - rj) is the 
dressed charge of the spin chain model [3l[ in zero 
external magnetic field. The parameter rj is related 
to the magnetic anisotropy of the exchange interaction 
cos r\= J z / J for easy-plane anisotropy, and cosh r\= J z / J 
for easy-axis anisotropy of the spin-spin interaction along 
the chain (A is 77-dependent). Hence, our experimental 
results permit a small anisotropy. However, an exact de- 
termination of 77 is limited by the experimental accuracy. 

In our experiments temperatures were small compared 
to the intrachain exchange interactions, T<CJ, J'. For low 
fields we can then estimate the temperature dependence 
of the relaxation rate A from the temperature dependence 
of the DSF by conformal field theory [3lj |: 



1 

A oc — — 

B V v 



(6) 



where a is the cut-off parameter, v=nJsm(r])/2r] is the 
spinon velocity v~irJ/2 for weak anisotropy, and 7 is 
related to the Luttinger liquid exponent, connected to 
the dressed charge. For the isotropic 5=1/2 AFHC we 
have 7«0. The temperature dependence of Ao in Fig. [2] 
is constant below approx. 150 K hence it is in agreement 
with 7«0. 

Formally the double chain in SrCu02 can be treated by 
introducing weak NNN exchange interaction. The main 
effect of the then frustrated NN and NNN exchange inter- 
action is an additional minimum (maximum) for q=n/2 
in the lower (upper) boundary of the DSF [32j. However, 
as we pointed out above, the main contributions to the 
/iSR (and NMR) spin-lattice relaxation rate have to come 
from q=0, 7T, which are not essentially changed by spin 



frustration in the chain. In addition, it is possible that 
the low temperature exponent 7 becomes non zero due to 
spin frustration t 33] . This effect should be small because 
we experimentally find 7«0 (see above). Furthermore, 
any effects of the NNN on the DSF would be detected by 
inelastic neutron scattering (INS). The absence of any 
such effect on the DSF measured by INS 1] shows that 
up to energies of the order of the exchange interaction 
J«2100 K the NNN interactions have negligible effects. 

It turns out that the results of Monte-Carlo simulations 
9] indicate diffusive behavior of the relaxation rate in the 
5=1/2 AFHC. Those results are valid for high and in- 
termediate temperatures i.e. T^>J, or TraJ. In our 
situation we have T«75 K, and J«2000 K, i.e. T<J, 
and we cannot apply the results of Grossjohann et al. [§] 
for the explanation of our experiments. Also, similar dif- 
fusive behavior was predicted within the field-theoretical 
approximation (lo| . However, the perturbation scheme 
that Sirker et al. use in Ref. [l(| for the calculation of 
the self-energy (mass operator) cannot be applied in our 
case, because both perturbation and the main part of the 
Hamiltonian are determined by the same exchange con- 
stant J, i.e. formally there is no small parameter in our 
situation. 

It is beyond the scope of this Letter to clarify the dif- 
ferences between our results and previous works 14- 
Here we only want to mention the importance of impu- 
rities. Hammerath et al. have recently shown 2l| that 
in SrCu02 bond disorder leads to the formation of a gap 
at low temperatures. This is not fully understood but 
it follows, that in experimental studies that probe the 
low energy spin excitations, impurities play a crucial role 
because they can alter the low energy excitation spec- 
trum away from that of the 5=1/2 AFHC. In most cases 
this circumstance is not as obvious because no gap is ob- 
served. However, it is clear that impurities can change 
the low energy excitation spectrum of the 5=1/2 AFHC 
drastically. The samples that we have studied are of high 
purity as has been shown by heat transport experiments 
[llj . They demonstrated that the mean free path of the 
magnetic heat conduction is of the order of 1 /jm 11 1. 
This is exceptionally high and can only be due to the 
high quality of our samples [ll| . We are therefore lead to 
believe that the apparent differences between our results 
and previous results 14 -111 are mostly due to the very 
low amount of impurities in our samples. 

In summary, the magnetic field dependence of the 
spin lattice relaxation rate A in SrCu02 is XocBq™ with 
n=0.9(3). This is in close agreement with n=l which fol- 
lows from the dynamic structure factor expected for the 
ground state of the isotropic antiferromagnetic 5=1/2 
Heisenberg chain (AFHC). We conclude that in this sys- 
tem the low energy spin dynamics is determined by the 
eigenstates of the Heisenberg spin chain Hamiltonian. 
Therefore, in SrCu02 spin transport is ballistic at low 
temperatures and fields 0<H, T-CJ. Furthermore, the 
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absence of a temperature dependence for T<150 K of the 
spin lattice relaxation rate A(0)=Ao in zero field indicates 
that frustration due to NNN interaction do not influence 
the low energy spin dynamics. In addition it shows that 
spin and heat transport [ll| are decoupled in this sys- 
tem. This docs not contradict the often cited magnetic 
contribution to the heat conductivity of low dimensional 
magnets. It follows, that the mechanism causing the tem- 
perature dependence of the the magnetic heat transport 
is not relevant for the spin transport in SrCu02- 
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